Search for: All records

Abstract Fault geometry affects the initiation, propagation, and cessation of earthquake rupture, as well as, potentially, the statistical behavior of earthquake sequences. We analyze 18,250 (−0.27 < M < 4.4) earthquakes of the 2016–2019 Cahuilla, California, swarm and, for the first time, use these high-resolution earthquake locations to map, in detail, the roughness across an active fault surface at depth. We find that the strike-slip fault is 50% rougher in the slip-perpendicular direction than parallel to slip. 3D mapping of fault roughness at seismogenic depths suggests that roughness varies by a factor of 8 for length scales of 1 km. We observe that the largest earthquake (M 4.4) occurred where there is significant fault complexity and the highest measured roughness. We also find that b-values are weakly positively correlated with fault roughness. Following the largest earthquake, we observe a distinct population of earthquakes with comparatively low b-values occurring in an area of high roughness within the rupture area of the M 4.4 earthquake. Finally, we measure roughness at multiple scales and find that the fault is self-affine with a Hurst exponent of 0.52, consistent with a Brownian surface. 

SUMMARY In the Gulf of California, Mexico, the relative motion across the North America–Pacific boundary is accommodated by a series of marine transform faults and spreading centres. About 40 M> 6 earthquakes have occurred in the region since 1960. On 2009 August 3, an Mw 6.9 earthquake occurred near Canal de Ballenas in the region. The earthquake was a strike-slip event with a shallow hypocentre that is likely close to the seafloor. In contrast to an adjacent M7 earthquake, this earthquake triggered a ground-motion-based earthquake early warning algorithm being tested in southern California (∼600 km away). This observation suggests that the abnormally large ground motions and dynamic strains observed for this earthquake relate to its rupture properties. To investigate this possibility, we image the rupture process and resolve the slip distribution of the event using a P-wave backprojection approach and a teleseismic, finite-fault inversion method. Results from these two independent analyses indicate a relatively simple, unilateral rupture propagation directed along-strike in the northward direction. However, the average rupture speed is estimated around 4 km s−1, suggesting a possible supershear rupture. The supershear speed is also supported by a Rayleigh wave Mach cone analysis, although uncertainties in local velocity structure preclude a definitive conclusion. The Canal de Ballenas earthquake dynamically triggered seismicity at multiple sites in California, with triggering response characteristics varying from location-to-location. For instance, some of the triggered earthquakes in California occurred up to 24 hr later, suggesting that nonlinear triggering mechanisms likely have modulated their occurrence. 

The interaction of aseismic and seismic slip before and after an earthquake is fundamental for both earthquake nucleation and postseismic stress relaxation. However, it can be difficult to determine where and when aseismic slip occurs within the seismogenic zone because geodetic techniques are limited to detecting moderate to large slip amplitudes or long duration small slip amplitudes. Here, we use repeating earthquakes (earthquakes that re‐rupture the same fault patch) as a proxy for aseismic slip during the 2011 Prague, Oklahoma earthquake sequence. We find that aseismic slip in the Prague earthquake sequence occurs both within the granitic basement and the overlying sedimentary rocks. The repeating earthquakes show that patches of aseismic slip are mostly located at fault intersections. These fault intersections hosted possible mainshock slip, abundant aftershocks, and afterslip. We estimate that ∼40% of the aftershocks are driven by afterslip. We interpret that aseismic slip occurs at fault intersections where stress heterogeneity creates patches of lower stress that are stable within a nonsteady state, rate‐state framework.

 

Megathrust earthquakes and their associated tsunamis cause some of the worst natural disasters. In addition to earthquakes, a wide range of slip behaviors are present at subduction zones, including slow earthquakes that span multiple orders of spatial and temporal scales. Understanding these events may shed light on the stress or strength conditions of the megathrust fault. Out of all types of slow earthquakes, very low frequency earthquakes (VLFEs) are most enigmatic because they are difficult to detect reliably, and the physical nature of VLFEs are poorly understood. Here we show three VLFEs in Cascadia that were dynamically triggered by a 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California. The VLFEs likely locate in between the seismogenic zone and the Cascadia episodic tremor and slip (ETS) zone, including one event with a moment magnitude of 5.7. This is the largest VLFE reported to date, causing clear geodetic signals. Our results show that the Cascadia megathrust fault might slip rapidly at some spots in this gap zone, and such a permissible slip behavior has direct seismic hazard implications for coastal communities and perhaps further inland. Further, the observed seismic sources may represent a new class of slip events, whose characteristics do not fit current understandings of slow or regular earthquakes.

 

Dynamic triggering of earthquakes has been reported at various fault systems. The triggered earthquakes are thought to be caused either directly by dynamic stress changes due to the passing seismic waves, or indirectly by other nonlinear processes that are initiated by the passing waves. Distinguishing these physical mechanisms is difficult because of the general lack of high‐resolution earthquake catalogs and robust means to quantitatively evaluate triggering responses, particularly, delayed responses. Here we use the high‐resolution Quake Template Matching catalog in Southern California to systematically evaluate teleseismic dynamic triggering patterns in the San Jacinto Fault Zone and the Salton Sea Geothermal Field from 2008 to 2017. We develop a new statistical approach to identify triggered cases, finding that approximately 1 out of every 5 globalM_w ≥ 6 earthquakes dynamically trigger microearthquakes in Southern California. The triggering responses include both instantaneous and delayed triggering, showing a highly heterogeneous pattern and indicating possible evolving triggering thresholds. We do not observe a clear peak ground velocity triggering threshold that can differentiate triggering earthquakes from nontriggering events, but there are subtle differences in the frequency content of the ground motion that may differentiate the earthquakes. In contrast to the depth distribution of background seismicity, the identified triggered earthquakes tend to concentrate at the edges of the seismogenic zones. Although instantaneously triggered earthquakes are likely a result of dynamic Coulomb stress changes, the cases of delayed‐dynamic triggering are best explained by nonlinear triggering processes, including cyclic material fatigue, accelerated transient creep, and stochastic frictional heterogeneities.

 

The vibrant evolutionary patterns made by earthquake swarms are incompatible with standard, effectively two-dimensional (2D) models for general fault architecture. We leverage advances in earthquake monitoring with a deep-learning algorithm to image a fault zone hosting a 4-year-long swarm in southern California. We infer that fluids are naturally injected into the fault zone from below and diffuse through strike-parallel channels while triggering earthquakes. A permeability barrier initially limits up-dip swarm migration but ultimately is circumvented. This enables fluid migration within a shallower section of the fault with fundamentally different mechanical properties. Our observations provide high-resolution constraints on the processes by which swarms initiate, grow, and arrest. These findings illustrate how swarm evolution is strongly controlled by 3D variations in fault architecture.